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The fractal structure of spin clusters and their boundaries in the critical two-dimensional Ising 
model is investigated numerically. The fractal dimensions of these geometrical objects are estimated 
by means of Monte Carlo simulations on relatively small lattices through standard finite-size scaling. 
The obtained results are in excellent agreement with theoretical predictions and partly provide 
significant improvements in precision over existing numerical estimates. 



I. INTRODUCTION 

The past few years have witnessed a surge in the ge- 
ometrical approach to phase transitions. The prototype 
of such an approach is percolation theory 0, which fo- 
cuses on clusters of (randomly) occupied sites or bonds 
on a lattice. The fractal structure of these geometrical 
objects and whether or not a cluster percolates the lattice 
are central topics addressed by the theory. Spin models 
such as the g-state Potts models can easily be mapped 
onto percolation theory, with neighboring spins in the 
same spin state lumped together in a cluster. Generally 
the geometrical spin clusters thus constructed do not per- 
colate at the critical temperature T c where the thermal 
phase transition takes place. But by erasing with a cer- 
tain temperature-dependent probability bonds between 
like spins, Fortuin and Kasteleyn (FK) showed that 
spin clusters can be constructed for the Potts models that 
do percolate at T c and encode the thermal critical behav- 
ior. They thus achieved a geometrical description of the 
thermal phase transition in these models. The cluster 
approach has been turned into an efficient Monte Carlo 
algorithm by Swendsen and Wang 3], and by Wolff 
where not individual spins are updated, as in local spin 
updates with, e.g., the Metropolis algorithm, but entire 
FK clusters. 

An exception to the rule that geometrical clusters do 
not percolate at T c is the two-dimensional (2D) g-state 
Potts model. The origin of this effect can be understood 
by extending the pure lattice model to include vacant 
sites. In a Kadanoff block-spin approach, such an exten- 
sion is natural as the vacant sites represent disordered 
blocks without a majority in any of the spin states, and 
is essential for establishing that the phase transition of 
the pure models changes from being continuous to first 
order at q = 4 0-0 In addition to the pure Potts critical 
behavior, the site diluted model also displays tricritical 
behavior at the same critical temperature T c . While the 
critical behavior of the pure model is encoded in the FK 
clusters, the tricritical behavior is encoded in the geomet- 
rical clusters 0H,@- Both cluster types percolate at T c . 
With increasing q, the critical and tricritical fixed points, 
which are characterized by the same central charge c, 
move together until merging at q = 4. 

Very recently, cluster boundaries of critical 2D systems 
have been studied analytically by means of the so-called 



stochastic Loewner evolution, introduced by Schramm 
[Tof . Various exact predictions for critical exponents pre- 
viously conjectured on the basis of the Coulomb-gas map 
|TTL ITa . fl3l j and conformal invariance 0] could rigor- 
ously be established by this method (for an overview see 
Ref. [3). 

In this paper, we numerically investigate the fractal 
structure of clusters in the 2D Ising model, correspond- 
ing to setting q = 2. We simulate the model on relatively 
small lattices (L = 8 — 512) with periodic boundary 
conditions, and apply standard finite-size scaling to de- 
termine the various fractal dimensions. In addition to 
studying the size or "mass" of FK and geometrical clus- 
ters, we also examine their boundaries. Those of geo- 
metrical clusters form the famous Peierls domain walls 
[l6|. separating spin clusters of opposite orientation. In 
a previous paper 0, we simulated these domain walls 
directly by considering the high-temperature represen- 
tation of the model. By duality, the high-temperature 
graphs, which are closed, are domain walls on the dual 
lattice. The closed graph configurations were generated 
by means of a Metropolis update algorithm, involving 
single plaquettes. 

Other recent numerical studies of the geometrical 
structure of 2D Potts models were reported in Refs. 0, 
ITU . Our results for the fractal dimensions are in excellent 
agreement with theoretical predictions [t! ll3L IT^. I20L l2l| . 
and provide in particular for the FK and geometrical clus- 
ters a considerable improvement in precision over the es- 
timates obtained in Refs. 0,0]- 

The rest of the paper is structured as follows. The next 
section summarizes the necessary theoretical background. 
Numerical results are presented in Sec. IIIII followed by 
concluding remarks in Sec. IIVI 

II. FRACTAL STRUCTURES 

The fractal properties of spin clusters and bound- 
aries, which are clusters themselves, are described by a 
straightforward extension of ordinary percolation theory 
0. Asymptotically, cluster distributions £ n take a gen- 
eral form 

£ n ~ n~ T exp(— On), (1) 

consisting of two factors: (i) an entropy factor, which 
measures the number of ways a cluster of size n can be 
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embedded in the lattice, and (ii) a Boltzmann weight, 
which suppresses large clusters when 9 is finite. Clus- 
ters proliferate and percolate the lattice when 9 tends 
to zero. The vanishing is characterized by an exponent 
a as 9 cx \T — T p | 1 /°', where T p denotes the percola- 
tion temperature. As explained in the Introduction, the 
percolation thresholds of both FK and geometrical clus- 
ters coincide with the thermal critical temperature of the 
2D Ising model. The entropy exponent r determines the 
fractal structure of the geometrical objects. Rather than 
extracting this exponent directly from the asymptotic 
behavior of a distribution at the percolation threshold, 
where the distribution becomes algebraic, it is expedient 
to extract it from derived quantities such as the perco- 
lation strength Poo, giving the fraction of sites in the 
largest cluster, and the average cluster size pj 
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(2) 



Since every site belongs to some geometrical and some 
FK cluster, the denominator in Eq. adds up to the 
total number of sites for these clusters. Close to the per- 
colation threshold, the observables obey the finite-size 
scaling laws j2^| 



(3) 



where L is the lattice size and £ the correlation length 
whose divergence at criticality is governed by the expo- 
nent v. The ratios /3/u and j/v are given in terms of r 
as 



i = d 3 - 



i' 



(4) 



with d the dimensionality of the lattice. The fractal di- 
mension D, which is also determined solely by the en- 
tropy exponent r, is related to these exponents via l] 



D 







T - 1 



(5) 



Generically, two (and only two) different cluster 
boundaries can be identified [23: the hull (H) and the 
external perimeter (EP), where the second can be un- 
derstood as a smoother version of the first. For 2D FK 
clusters, the two boundaries are in one-to-one correspon- 
dence, with their fractal dimensions satisfying the rela- 
tion pn 



(6) 



The map transforming one FK boundary dimension into 
the other conserves the central charge c, which may be 
parametrized as [I^.IT3| 



, 6(1 - S) 2 „ „ / 1 
c = 1 - v _ ' = 13 - 6 I k + - 

K \ AC 



(7) 



where 2 > k > 1 parametrizes the two-dimensional q- 
state Potts models 



y/q = —2 cos(7t/k), 



(8) 



with < q < 4. In terms of k, the fractal dimensions of 
the FK boundaries 0, 0] can be expressed as 



n FK 



n FK 



1 



(9) 



while the central charge conserving map corresponds to 
letting k — > 1/k. These explicit forms are seen to sat- 
isfy the duality relation 10. With the scaling relations 
(J5J, the critical exponent ratios characterizing the FK 
boundaries become 



7^ ^ 



7lp> = 1/R, 



(10) 



where a single correlation length with exponent v is as- 
sumed. It thus follows that the two FK boundary sizes 
scale with inverse exponents: 



v FK 

Xh 



v FK 

Xep 



(11) 



In contrast to FK clusters, geometrical clusters are 
characterized by only one boundary dimension, i.e., the 



fractal dimensions £>§ and D^. p of the hull and external 
perimeter coincide, D§ 



Such cases ar e sig naled 
by a negative fractal dimension of the red sites [I3 , sites 
that, when removed, lead to a splitting of the cluster into 
disconnected parts. 

The central charge conserving map R — > 1/k trans- 
forming the hull dimension of FK clusters into that 
of their external perimeters, -Dgp, also maps it onto the 
hull dimension of geometrical clusters, implying pi) 



LJ Ep — U H . 



(12) 



This relation is remarkable as it involves the two different 
boundary types. In the context of uncorrelated percola- 
tion j23|, the hull of a cluster in 2D is traced out by a 
directed random walker constrained to move on the clus- 
ter only, whereas the external perimeter is traced out by 
a walker constrained to move around the hull on sites 
neighboring, but not belonging to, that cluster. 

To numerically verify relation (|12(l . we wish to treat 
external perimeters of FK clusters and hulls of geomet- 
rical clusters in the same manner. To this end we apply 
the same algorithm used to trace out geometrical hulls to 
find the FK external perimeters, where it is recalled that 
FK clusters differ from geometrical clusters in that with a 
prescribed temperature-dependent probability bonds are 
erased. The difference between tracing out FK hulls and 
external perimeters then reduces to (see Fig. ^) allowing 
the random walker to move to a nearest neighbor site on 
the FK boundary only provided the connecting bond is 
set (hull) or always (external perimeter). 

To conclude this section we list in Table|Jthe predicted 
exact values 0, 0, 0, H3, OHJ for the various fractal 
dimensions and corresponding critical exponents we wish 
to determine numerically. 
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FIG. 1: (Color online) Part of a single FK cluster of nearest 
neighbor sites (filled circles) connected by bonds (black links) . 
(a) Sites belonging to the hull (dark filled circles) are found 
by allowing the random walker tracing out the boundary to 
move only over set bonds, (b) Sites belonging to the exter- 
nal perimeter (dark filled circles) are found by allowing the 
random walker to move to a nearest neighbor on the cluster 
irrespective of whether the connecting bond is set or not. The 
external perimeter, which contains two sites less than the hull 
for this boundary segment, is therefore a smoother version of 
the hull. 



ters, as well as of the FK hulls and external perimeters. 
Where j/u and j3/v are listed, both are measured in- 
dependently by considering the average cluster size \i 
which gives ^jv according to Eq. © with X(0) = const, 
and the percolation strength P oai which gives f3/u accord- 
ing to Eq. © with P(0) = const. The data were fitted 
using the least-squares Marquardt-Levenberg algorithm. 

While including the percolating clusters when consid- 
ering the mass of the cluster, we ignore them in tracing 
out cluster boundaries. Because of the finite lattice size, 
large percolating clusters have abnormal small (external) 
boundaries, so that including them would lead to a dis- 
tortion of the hull distribution. Moreover, the Grossman- 
Aharony algorithm |23| we use to trace out the cluster 
boundaries generally fails on a percolating cluster as its 
boundary not necessarily forms a single closed loop any 
longer. With the percolating clusters ignored, the bound- 
ary exponents /3h,ep/^ cannot be determined, while the 
summation in the expression J2J) for the average cluster 
boundary size is restricted to nonpercolating clusters. 







7h/^ 


-Dep 7ep/^ 


FK 


15/8 7/4 1/8 5/3 


4/3 


11/8 3/4 


G 


187/96 91/48 5/96 11/8 


3/4 





TABLE I: Predicted values for the fractal dimensions with 
the corresponding critical exponents characterizing Fortuin- 
Kasteleyn (FK) and geometrical (G) clusters (C), their hulls 
(H), and external perimeters (EP). 



III. MONTE CARLO SIMULATIONS 

The simulation data was collected on square lattices 
of linear size L = 8, 10, 12, 14, 16, 20, 24, 32, 40, 48, 64, 90, 
128, 180, 256, 360, and 512 with periodic boundary condi- 
tions, using the Swendsen-Wang cluster algorithm [3j in 
about 5 x 10 4 measurements at the critical temperature 
T c = 2/ln(l + V2), every rth sweep of the lattice, where 
t denotes the autocorrelation time (rounded off to the 
next largest integer). The value of r for the various lat- 
tice sizes L was estimated from the energy time series to 
vary from r sw 4 for L = 8 to r w 9 for L = 512. We have 
chosen the energy time series here as it generally leads 
to a conservative estimate of the autocorrelation time for 
cluster algorithms. A total of 5 x 10 3 lattice sweeps were 
used for equilibration. Statistical errors were estimated 
by means of jackknife binning |24| . 



A. Fractal Dimensions 

Tables HTl 11111 II V I andlVlsummarize the values obtained 
for the critical exponents of the FK and geometrical clus- 



1. FK Clusters 

In Ref. HI. the value Dg K = 1.87(1) compared to the 
exact result [19] D FK — 15/8 = 1.875 was reported for 
the fractal dimension of FK clusters. It was obtained on a 
single, but very large lattice (L — 2 12 = 4096) with both 
open and periodic boundary conditions by considering 
the number of bond clusters as a function of the radius 
of gyration. The authors observed a slow and complex 
approach to the asymptotic behavior and therefore in- 
cluded corrections to scaling to arrive at their numerical 
estimates. To limit the number of fit parameters, they 
scanned the fractal dimensions in the neighborhood of 
the predicted values, and completely fixed the values of 
the correction exponents to the theoretically predicted 
ones. This left them with still four parameters to fit. 
Error bars on the values of the fractal dimensions were 
determined as the range where x 2 /d.o.f. < 2. 

In this study, where we use different lattice sizes and 
consider not the bonds, but the sites in a cluster, we find 
a simple approach to the asymptotic behavior with very 
small corrections to scaling in the observables required to 
determine D^ K all the way down to the smallest lattice 
considered (L = 8). We can therefore apply finite-size 
scaling without correction terms. Since the fits involve 
only two parameters, no exponents need to be fixed be- 
forehand. To minimize the effect of unavoidable (small) 
corrections to scaling, we pick the fit over the largest lat- 
tice sizes given in Table HT1 i.e., over the range 64 — 512, 
leading to 

D F C K = 1.8753(6) (13) 

with x 2 /d.o.f = 0.61 from the average cluster-size data 
and 

D F C K = 1.8752(8) (14) 



4 



Fit Interval 


7c7" 


X 2 /d.o.f. 




X 2 /d.o.f. 


8 - 256 


1.7512(6) 


0.89 


0.1244(4) 


0.74 


16 - 256 


1.7507(8) 


1.02 


0.1246(5) 


0.95 


32 - 256 


1.7500(12) 


1.39 


0.1249(8) 


1.31 


40 - 256 


1.7507(13) 


1.43 


0.1246(9) 


1.43 


64 - 256 


1.7496(21) 


0.52 


0.1254(14) 


0.33 


8 — 512 


l.75ll(4) 


0.86 


0.1244(3) 


0.71 


to — oiz 


1. (out (0) 


U.90 


n 1 r iAafA\ 
U.lz40^4J 


A on 

U.oB 


32 - 512 


1.7505(8) 


1.23 


0.1247(5) 


1.14 


40 - 512 


1.7509(9) 


1.20 


0.1245(6) 


1.17 


64 - 512 


1.7505(12) 


0.61 


0.1248(8) 


0.47 


exact 


1.7500 = 7/4 




0.1250 = 1/8 





TABLE II: FK clusters. 



with x 2 /d.o.f = 0.47 from the percolation-strength data. 
Both estimates are well within one standard deviation 
from the exact prediction [ig| £>£ K = 15/8 = 1.875. 

A possible explanation for the improved accuracy we 
achieved over Ref. [l8], although working on smaller lat- 
tices, may be that in that study clusters touching the 
boundary were ignored in the measurements. This set in- 
cludes all the percolating clusters. As will be illustrated 
in the next subsection, omitting percolating clusters can 
lead to strong corrections to scaling. 



2. Geometrical Clusters 

In Ref. [13, the values -yg/v = 1.901(11) and /3g/zy = 
0.052(2) were reported for geometrical clusters. These 
results were obtained on lattices ranging in size from 
L = 600 to 2000, i.e., again much larger than the ones 
considered by us. Instead of using periodic boundary 
conditions, as we did, free boundary conditions were 
adopted. Another difference from our ap proa ch is that 
percolating clusters were excluded in Ref. \u\ . 

Our estimates, obtained from the largest lattice sizes 
listed in Table ITTT1 i.e., from the interval 64 — 512, are 



U Q 



1.9476(3) 



(15) 



with x 2 /d.o.f = 0.44 from the average cluster-size data 
and 



D% = 1.9473(4) 



(16) 



with x 2 /d.o.f = 0.23 from the percolation-strength data. 
Both are in excellent agreement with the exact prediction 
D% = 187/96 = 1.9479 . . .. 

In Fig. |21 we show our data for the average cluster size 
obtained by including all clusters and compare it with the 
data obtained with percolating clusters excluded from the 
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FIG. 2: (Color online) Log-log plot of the average size Xc °f 
critical geometrical clusters as a function of the lattice size 
L. Statistical error bars are smaller than the symbol size in 
the figure. The straight lines proportional to L 91 / 48 are put 
through the data points at L = 512 by hand to demonstrate 
the corrections to scaling for smaller lattice sizes when per- 
colating clusters are excluded, and the absence thereof when 
they are included in the measurements. 



Fit Interval 


7c lv 


X 2 /d.o.f. 




X 2 /d.o.f. 


8-256 


1.8941(3) 


1.16 


0.0532(2) 


0.83 


16 - 256 


1.8944(4) 


1.39 


0.0530(3) 


0.94 


32 - 256 


1.8946(5) 


1.90 


0.0528(4) 


1.24 


40 - 256 


1.8950(6) 


1.78 


0.0526(4) 


1.09 


64 - 256 


1.8949(9) 


0.47 


0.0527(7) 


0.18 


8-512 


1.8943(2) 


1.29 


0.0531(2) 


0.89 


16 - 512 


1.8946(3) 


1.34 


0.0529(2) 


0.88 


32 - 512 


1.8949(4) 


1.60 


0.0528(3) 


1.02 


40 - 512 


1.8951(4) 


1.39 


0.0526(3) 


0.86 


64 - 512 


1.8951(5) 


0.44 


0.0527(4) 


0.23 


exact 


1.8958 = 91/48 




0.0521 = 5/96 





TABLE III: Geometrical clusters. 

measurements, similarly to what was done in Ref. 
While virtually absent in the former, corrections to scal- 
ing are present in the latter case. This may explain why, 
although working on smaller lattices, we obtained much 
better estimates than in Ref. [17|. We found similar cor- 
rections to scaling when instead of percolating clusters, 
the largest cluster in each measurement was excluded, as 
is common in random percolation . 



3. FK Hulls 

Surprisingly, the results for the hulls of FK clusters 
jiven in Table ITVl show a clear tendency to the predicted 
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FIG. 3: (Color online) Log-log plot of the average hull size 
Xh K of critical FK clusters as a function of the lattice size L. 
Statistical error bars are smaller than the symbol size in the 



figure. The straight line 0.310 L 
two-parameter fit in the interval 8 



is obtained through a 



48. 



value [13j 7h K / j/ = 4/3, corresponding to = = 
1.6666 . . ., when restricting the fitting window to increas- 
ingly smaller lattice sizes. For example, for the interval 
8 - 48 we hnd 



1.665(3) 



(17) 



with x 2 /d.o.f = 0.79, indicating a good fit. This es- 
timate, which should be compared with the estimate 
1.66(1) given in Ref. [l^, is within one standard devi- 
ation from the exact prediction. From Fig. [3]we see that 
the average FK hull size measured on larger lattices falls 
below the expected value extrapolated from smaller lat- 
tices. We have not been able to determine the cause of 
this behavior. In fact, when fitting not at the low end 
but at the high end of the lattice sizes considered, we 
obtain fits of comparable quality, but the estimate for 
the exponent converges to a value significantly below the 
predicted one (see top part of Table II V B . 



4- Geometrical Hulls 

As for the clusters' mass when disregarding percolat- 
ing clusters, we observe strong corrections to scaling for 
the hulls of geometrical clusters (see Fig.0J. This is dif- 
ferent from what we found using the plaquette update 
to directly simulate the hulls of the spins on the dual 
lattice 9J, where these corrections were virtually absent 
(see Fig. 11 of that paper), allowing us to obtain precise 
estimates for the critical exponents on relatively small 
lattices. In that study, the largest hull was omitted in 
each measurement. 

To understand the strong corrections to scaling found 
here for the geometrical hulls, we depict in Fig. [5] the 
geometrical cluster and the corresponding hull distribu- 
tions for L — 32 normalized to the volume L 2 . The bump 
at the tail of the cluster distribution is due to the finite 
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1.281(5) 


0.62 
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1.311(2) 


2.88 


8- 


128 


1.318(3) 


1.50 


8- 


90 


1.321(3) 


1.29 


8- 


64 


1.325(4) 


1.10 


8- 


48 


1.329(5) 


0.79 



exact 



1.333 = 4/3 



TABLE IV: Hulls of FK clusters. 



Fit Interval 


7IfV" 


X 2 /d.o.f. 


24- 


256 


0.763(3) 


4.32 


32- 


256 


0.758(3) 


3.10 


40- 


256 


0.755(4) 


3.24 


64- 


256 


0.748(6) 


4.13 


24- 


512 


0.752(2) 


7.97 


32- 


512 


0.747(2) 


5.73 


40- 


512 


0.744(3) 


5.29 


64- 


512 


0.736(4) 


4.37 



exact 



0.750 = 3/4 



TABLE V: External perimeters of FK clusters. 



size of the lattice, with percolating clusters gulping up 
smaller ones reached by crossing lattice boundaries. The 
subsequent sharp drop-off arises because of the limited 
number of lattice sites available. 

Initially, as Fig. [S] clearly shows, the hull distribution 
follows more or less the cluster distribution. This is a 
common feature of all boundary distributions considered. 
The relatively early drop-off of the hull distribution is 
because we omit percolating clusters when tracing out 
cluster boundaries. As a result, the average hull size is 
underestimated and the data points in Fig. 01 are below 
the expected line extrapolated from larger lattice sizes. 
With increasing lattice size, the effect becomes smaller 
(see Fig. where also the distributions for L = 512 
are included) and the data points approach the expected 
asymptotic scaling, corresponding to [2(j 7h/^ = 3/4, 
and L>g = 11/8. Fi gure [o] shows in addition that the 
asymptotic behavior of the hull distribution sets in for 
relatively large hull sizes (n > 100). On smaller lat- 
tices, the asymptotic behavior can therefore simply not 
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FIG. 4: (Color online) Log-log plot of the average hull size 
X§ °f critical geometrical clusters as a function of the lattice 
size L. Statistical error bars are smaller than the symbol 
size in the figure. The straight line proportional L 3 ^ 4 is put 
through the data point at L — 512 by hand to demonstrate 
the strong corrections to scaling for smaller lattice sizes. A 
three-parameter fit in the interval 16 — 512 gives the value 
w = 0.54(3) for the effective correction-to-scaling exponent. 



be probed, explaining the strong corrections in Fig. 

To see if our data are at least consistent with the the- 
oretical prediction, we account for corrections to scaling 
by fitting the average hull-size data to the form 



X % = aL^(l-bL-™) 



(18) 



with an effective correction-to-scaling exponent w. We 
fix 7§/^ = 3/4 to the predicted value, leaving us with 
three parameters (w and the two amplitudes a and b) to 
fit. For the interval 16 — 512 we obtain 



= 0.54(3) 



(19) 



with % 2 /d.o.f = 1.60, indicating a reasonable fit (see 
Fig.^J and therefore consistency with the theoretical pre- 
diction for D^. 



FIG. 5: (Color online) Log-log plot of the (normalized) ge- 
ometrical cluster and hull distributions at T c on the largest 
lattice considered (L = 512) and on a relatively small lattice 
(L = 32). The number of measurements taken on the largest 
lattice was about 5 x 10 4 as in most part of this paper. Sta- 
tistical error bars are omitted from the data points for clarity. 
On the smaller lattice, about 5 x 10 5 measurements, which 
is an order of magnitude more than used in the rest of the 
paper, were taken to achieve good statistics. Here, the statis- 
tical error bars are smaller than the symbol size in the figure. 
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5. FK External Perimeters 

The corrections to scaling are less pronounced for the 
external perimeters of FK clusters as they are generally 
smaller than geometrical clusters and thus less likely to 
percolate (see Fig. EJ. Also the asymptotic behavior is 
reached earlier than for geometrical hulls. The smallness 
of the corrections allows us to obtain reasonable fits for 
the average external perimeter size. Our result obtained 
from the fitting interval 64 — 512 in Table ]V\ yields the 
fractal dimension 

£>|£ = 1.368(2) (20) 

with x 2 /d.o.f — 4.37. This estimate is compatible with 
the exact prediction [2lJ = 11/8 = 1.375, and im- 

proves the estimate D^p = 1.36(2) reported in Ref. 
by about one order of magnitude. 



FIG. 6: (Color online) Log-log plot of the average external 
perimeter size Xep °f FK clusters at T c as a function of the 
lattice size L. Statistical error bars are smaller than the sym- 
bol size in the figure. The straight line 1.388 L ' 736 is obtained 
through a two-parameter fit in the interval 64 — 512. Note the 
corrections to scaling for smaller lattice sizes. 



B. Distributions 

In Fig. [?1 the distributions of the three different bound- 
aries studied are plotted for L — 512 to show the 
crossover of the external perimeters of FK clusters. Start- 
ing similarly to the FK hull distribution, the FK external 
perimeter distribution asymptotically approaches that of 
the geometrical hulls, in accordance with relation (|12H . 
In other words, the FK external perimeter distribution 
interpolates between that of the FK (for small n) and 
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FIG. 7: (Color online) Log-log plot of the distributions of 
the three different boundaries at T c for L — 512. Statistical 
error bars are omitted from the data points for clarity. The 
FK external perimeter distribution initially follows the FK 
hull distribution before at around n ~ 100 it crosses over to 
its asymptotic behavior which it shares with the geometrical 
hulls. 

geometrical hulls (asymptotically for large n). 

Figure [S] summarizes all the cluster and boundary dis- 
tributions studied for L = 16 and 512. The distributions 
are normalized to the volume L 2 . Upon increasing the 
lattice size, the normalized distributions tend to a univer- 
sal curve. The slow approach to the asymptotic form of 
the geometrical hull distribution (and to a lesser extent 
that of the FK external perimeter distribution), with the 
associated strong corrections to scaling we observed for 
these objects, stands out clearly from the other distribu- 
tions. 

It is tempting to directly analyze the distributions £ n 
and to extract the exponent r from the asymptotic be- 
havior of £ n , which is algebraic at criticality. However, 
this method gives far less accurate results than apply- 
ing finite-size scaling to observables involving the sum 
over the cluster sizes n. The main drawback of the 
method is the great sensitivity to the location of the fit- 
ting window, i.e., the interval of n. The fitting range 
cannot be started at too small cluster sizes, where the 
distribution has not taken on its asymptotic form yet, 
while too large cluster sizes, which are generated only a 
few times during a complete Monte Carlo run, are also 
to be excluded because of the noise in the data and finite 
size effects. Even in those cases for which we obtained 
very accurate results through finite-size scaling analyses, 
the direct analysis of the distributions gave unsatisfac- 
tory results. 

IV. CONCLUSIONS 

In this paper, the fractal dimensions of spin clusters 
and their boundaries appearing in the critical 2D Ising 



model were studied numerically. The Monte Carlo simu- 
lations were carried out on comparatively small lattices. 
Standard finite-size scaling was applied to obtain very 
precise estimates for the cluster dimensions, significantly 
improving existing ones. The results confirm the exact 
theoretical predictions to a high degree of precision. 

For the boundary dimensions, although improving ex- 
isting estimates, we obtained less accurate results because 
of corrections to scaling. We observed the strongest cor- 
rections for the geometrical hulls, whose distribution ap- 
proaches its asymptotic form very slowly. In a previous 
numerical investigation j9( , where we simulated the hulls 
of the spins on the dual lattice directly, corrections to 
scaling were virtually absent, allowing us to establish the 
geometrical hull dimension to fairly high precision. 

To our surprise, we found the fractal dimension of the 
FK hulls to converge to the predicted value only when 
restricting the fitting window to increasingly smaller lat- 
tice sizes. In general, one expects of course such a con- 
vergence when increasing the lattice size, rather than de- 
creasing it, so as to minimize corrections to scaling. The 
measured average FK hull size on larger lattices falls be- 
low the line extrapolated from smaller lattices. The cause 
for this behavior eludes us. 

To verify relation (|12fl . involving the two different 
boundary types that can be defined for a cluster, viz. 
hulls and external perimeters, we treated the two bound- 
ary types in a similar manner. Usually, hulls are traced 
out by a directed random walker on the cluster whereas 
external perimeters are traced out by a directed random 
walker probing the cluster from the outside. We, on the 
other hand, applied the hull algorithm also to the exter- 
nal perimeters of FK clusters, with the proviso that the 
random walker can move to a nearest neighbor site on 
the FK boundary even when the connecting bond is not 
set (for the hull, the bond must be set). 
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FIG. 8: (Color online) Log-log plot of the critical distributions studied for L — 16 (short curves) and 512 (long curves), (a) 
The (normalized) distributions of FK clusters and of their hulls and external perimeters. For clarity, the latter two are shifted 
downward by two decades each, (b) The (normalized) distributions of geometrical clusters and of their hulls. The latter is 
shifted downward by four decades. Statistical error bars are also for clarity omitted from the data points. The straight lines 
are obtained through one-parameter fits with the slopes fixed to the predicted values. To achieve good statistics for the L — 16 
lattice, about 5 x 10 5 measurements were taken-an order of magnitude more than used for the L = 512 lattice. 
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